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Abstract 

This paper aims to introduce a new parameterisation for the coupling Q in interacting dark matter and 
dark energy models by connecting said models with the Continuous Tower of Scalar Fields model. Based 
upon the existence of a dark matter and a dark energy sectors in the Continuous Tower of Scalar Fields, a 
simplification is considered for the evolution of a single scalar field from the tower, validated in this paper. 
This allows for the results obtained with the Continuous Tower of Scalar Fields model to match those of an 
interacting dark matter - dark energy system, considering that the energy transferred from one fluid to the 
other is given by the energy of the scalar fields that start oscillating at a given time, rather than considering 
that the energy transference depends on properties of the whole fluids that are interacting. 


1. Introduction 

One of the deepest problems in Cosmology lies in the lack of understanding of about 95% of the energy 
content of the universe, the so called dark sector [l|. 

Composing this sector, two components with clearly distinct behaviours exist, dark matter (DM) and 
dark energy (DE), with the energy density of DM being around half the energy density of DE today. 
While DM behaves like matter, lacking only the capability to interact electromagnetically, DE behaves 
like a cosmological constant, a fluid with constant energy density (details on why considering a simple 
cosmological constant would present some problems can be found in [2|). 

Scalar fields are generally good candidates for both components of the dark sector, as it is possible to 
replicate the same behaviour of DM or DE when an adequate potential for the scalar field is chosen. While 
most DE models consider scalar fields Q, research on DM spans from Cosmology to Particle Physics. On 
the cosmological side, there are several models considering from Axions to discrete towers of scalar fields 
0, with the potential considered for the scalar fields in both mentioned cases being V = Despite 

most models consider scalar fields to explain a single component of the dark sector, there are some models 
that aim to account for both dark components, like the Continuous Tower of Scalar Fields (CTSF) model 
presented in Q. 

Another field of study that is particularly interesting when it comes to the dark sector lies in possible 
interactions between DM and DE, which is an obvious possibility when taking into account that the densities 
of both components are quite similar today, as mentioned before. 

The fact that those energy densities are of the same order is the so-called ” coincidence problem”, despite 
many argue it is not really a problem (a good discussion on this topic can be found in [7[). This possible 
interaction between the two components of the dark sector could possibly alleviate this problem, and explain 
why the energy densities of both components are quite similar today. 
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The most standard way to model this interaction consists in considering both DM and DE as fluids and 
add an interaction term to the contiiruity equations of both fluids in the following way: 


Pde + (1 + w)pde = Q 

Pdm + ^HpDM = —Q , 


( 1 ) 


The first parameterisations for the coupliirg term Q were given by Q = 5dmH poM, [3; which was later 
generalised to Q = SdeHpde + SdmH poM, l9|. In said parameterisations, the 6 terms give the strength 
of the coupling, which can either be constant or variable in time. By taking into account that the ratio 
R = Pde/Pdm today must agree with the observations, it is possible to constraint 5, which must obey the 
followirrg relation for the models mentioned, |(5| << 1. Also, it is relevant to mention that, rather than 
considerirrg a 5H factor, some models chose to describe the strength of the interaction by a E factor instead. 

With this formalism, it is necessary to resort to dynamical variables and a dynamical systems kind of 
analysis in order to understand how the energy densities of DE and DM evolve in time, given that it is not 
possible to solve the differential equations presented in [T] independently. 

It is also relevant to mention that DM and DE do not necessarily need to be fluids, and in fact is also 
possible to consider the interaction between a fluid and a scalar held, the so-called interacting quintessence 
model [13, or even the interaction between two scalar Helds [llj |. in which the interaction is achieved by 
considering a modihed potential that contains the coupling between the two scalar Helds. 


In this paper, the parameterisation for the coupling Q considered will be considerably different to the 
ones presented before. 

Rather than looking at DM and DE as just Huids, those two components will be also looked as sets of 
scalar Helds that mutate from a DE state into a DM state, and hence the parameterisation of Q will describe 
this transference of scalar Helds. Rather than depending on the total energy of the DM and DE fluids, as 
done before, Q will depend on the energy of the scalar Helds that shift from one state into the other. 

The physical motivation for such a parameterisation for Q lies in the CTSE model, in which scalar Helds 
start by behaving as DE, with said behaviour shifting into DM at later times. 

Therefore, the paper will develop as follows. In section 2, a simpliHed CTSE model, considering a 
simplifled behaviour for the single scalar Helds that compose the CTSE is presented and validated, while in 
section 3, the interacting DM-DE equations are used in order to obtain the same results for the evolution of 
the CTSE as in section 2. The paper ends with section 4, which features the flnal discussions of the results 
obtained in the paper. 


2. Simplified Continuous To-wer of Scalar Fields 

As mentioned in the introduction, scalar Helds are usually one of the best candidates to explain the dark 
sector, given that they can behave as DM or DE, depending on the potential chosen. 

In (6|, a new innovative model was proposed to unify the dark sector. To obtain this uniflcation, a set 
of minimally coupled scalar Helds with V = and a continuous distribution of masses is considered. 

Since no coupling is considered between the scalar Helds in said set, each scalar field will independently obey 
to the usual Klein-Gordon equation: 


+ = 0 ( 2 ) 

and the energy density for each scalar field will therefore be given by: 

P'Pm. = 


2 


( 3 ) 



With the potential mentioned before, the Klein-Gordon equation presented in eqU] is equivalent to the 
equation to a damped harmonic oscillator, with wq = m. 

Considering the previous analogy, it becomes clear that it is possible to choose a set of initial conditions 
for which the scalar fields are constant until t Ki T, oscillating afterwards. 

Therefore, the scalar field will behave like DE until it starts oscillating. Afterwards, the behaviour will 
depend on how H evolves in time as well, and it will therefore depend on the era considered. 

During radiation and matter domination, H can be parameterised as H = p/t, with p = 1/2 and p = 2/3 
respectively, which means the energy of the scalar field will then decay with a~^ oc behaving like DM. 

When the CTSF dominates, the situation becomes much more complex. Generally, it is not possible to 
parameterise H with a simple expression before solving the Klein-Gordon equation, due to the fact that 
H will now depend on the energy densities of the individual scalar fields on the CTSF, which means those 
scalar will in fact be now coupled gravitationally. This leads to two distinct cases for the behaviour of the 
individual scalar fields, which will depend on which kind of scalar fields that dominates. 

If the CTSF is dominated by the scalar fields that are already oscillating, then H can be parameterised by 
H = p/t, with p = 2/3, just like during matter domination, and hence the evolution of the individual scalar 
fields will be the same as during said era, behaving like DM. This case is equivalent to matter domination 
with DM being given by the sector of the CTSF that contains the oscillating scalar fields, and therefore the 
results presented for matter domination will be valid in this case as well. 

On the other hand, if most of the energy of the CTSF is due to the scalar fields that are frozen, the 
behaviour of the individual scalar fields will be considerably different. When t k, T, the scalar fields will 
start experiencing over-damped oscillations for which k, 0. This means that, despite the scalar fields are 
oscillating, their behaviour is actually closer to DE than to DM, despite it is still something in between those 
two regimes. The longer this regime lasts, the longer the scalar fields experience over-damped oscillations, 
which means the behaviour of a single scalar field is hard to parameterise, but also that the transition from 
a DE state to a DM state cannot be approximated by an instantaneous transition, as done in the previous 
eras. 

It is also relevant to mention that the asymptotic behaviour of the scalar fields is a DM behaviour, 
which means that if a given scalar field is already behaving like DM, its behaviour will not change when the 
universe enters an era dominated by the CTSF. 

For the two simple cases mentioned before, radiation and matter domination domination, the total 
evolution of can therefore be approximately described by two separate cases. If the scalar field starts 
oscillating during radiation domination, its evolution can be roughly given by: 


P<^7n — P'Pm^i'ni) ^ ^ 


1 , t <td 

, td <t <te 

, te < t 


(4) 


with td = l/m, which approximates the time for which the scalar fields start oscillating, given that cjq = m 
and ff is of order 1/t, and tmi, which is the initial time, while t^ is the time of the transition from radiation 
to matter domination, which means it will be relevant for m > l/tg. If the scalar field has m < 1/te, it will 
not start oscillating before matter domination and so its evolution can be approximately described by: 


Pfp-m — P4>m(/''^ni) ^ 


1 


, t <td 
, td<t 


(5) 


Physically, this can be seen as if the scalar fields have only the possibility to be in two different states. 
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a DE one (when the scalar fields are frozen), and a DM one (when the scalar fields oscillate), with the 
transition from one state to the other being instantaneous at t = td- 

However, due to the reasons mentioned before, this kind of simplification is only valid until the DE sector 
of the CTSF starts to dominate. 

For the evolutions of the scalar fields, described by eqsll]and[5l to happen, the following initial conditions 
must be considered: 


4^7n(j'ini') — | m I 

4^m — 0 5 


( 6 ) 


as presented in 


With said initial conditions, the initial energy density for the scalar field is given by: 


(7) 

with C = since the potential term in the energy dominates. 

In the previous expressions, A is the amplitude for the scalar field, cr is the width of the distribution for 
the CTSF (which can be roughly translated into the highest mass considered in the CTSF), m the mass of 
a scalar field in the CTSF and b a parameter that allows the energy density to be shaped over the mass 
(6=1 allows for the energy of the scalar fields to be independent of the mass while for 6 < 1 the energy 
grows with the mass and for 6 > 1 the opposite occurs). 

It is relevant to mention that the previously presented expressions for the initial conditions, originally 
presented in Q, are quite general given said paper aimed to study the implications of the CTSF without 
paying special attention to the underlying mechanism generating it. This resulted in a more general model, 
but it is still a model with proper physical motivation, given that the CTSF can in fact be a good approxi¬ 
mation to a discrete tower with a large number of scalar fields and a small mass gap. Therefore, unparticle 
physics [l^ , braneworld scenarios (like the Randall-Sundrum [l^ or DGP [l^ models), or even more exotic 
setups could provide the underlying mechanism to generate the CTSF (a deeper discussion on the subject 
is presented in the introduction of [6|). 


By plottiM the simplified solutions, presented in eqsHl and [SJ together with the analytic solutions, 
obtained in (6|, it is possible to notice that the only significant differences in behaviour occur around 
t = 1/to, as expected. However, afterwards the value of is slightly shifted during radiation domination, 
which does not occur for matter domination, as shown in figure [TJ 

Given said solutions provide good approximations to the analytic solutions, it is now possible to use 
them to obtain a new simplified approximation solutions for the energy density of the whole CTSF, px, 
during radiation and matter domination. 

However, it is relevant to point out another simplified solution for px, during radiation and matter 
domination, was already presented in [^, but as will be seen in the next section, the new set of solutions 
obtained afterwards will be crucial to show the connection between the CTSF and models with interacting 
DM-DE, as this paper aims to. 

The process to determine this new simplified solution for px is similar to what was done in Q, but to 
keep the calculations very simple, further simplifications to the CTSF model were considered. 

As mentioned in Q , an Heaviside step function can be a good approximation to the Gaussian distribution 
and will therefore be used in this study. The only parameter of the distribution will therefore be cr, which 
is the higher mass considered in the distribution, so that all scalar fields with to < cr are included in the set 
with the same weight (with the energy density of each scalar field depending only on the initial conditions). 

With the set of scalar fields defined, it must now be taken into account that there are two different 
sectors present in the GTSF, one in which the scalar fields are frozen and the other in which the scalar fields 
are oscillating, with the mass that separates both sectors being given by to^ = 1/t, which will be simpler 
than considering rus = H, as done in [6|, for the aimed calculations. 
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Figure 1; In the left panel, is plotted during radiation domination, while the right plot shows the same quantity during 
matter domination. In both plots, the thick line denotes the approximation considered in this section, while the dashed lined 
refers to the analytic solution obtained in In both cases, the energy density decays like matter, with a“®. 


Given that the frozen scalar fields will behave as DE, while the oscillating scalar fields will behave as 
DM, pt for the CTSF can be obtained by summing the energy for the two mentioned sectors: 

Pt = Pdm + Pde (8) 

As seen before, pde does not depend on whether the CTSF is in radiation or matter domination. 
Assuming that tini = l/cr, pde is given by the following integration of over the scalar fields: 

nTTis ry 

PDE = ^* J ^ 

with G* = Cja = in order to normalize the energy density. The energy lost by this sector of the 

CTSF is only due to the most massive scalar fields start oscillating, when t = 1/m, and therefore passing 
to the DM sector. 

Unlike the evolution of pde, the evolution of pdm will depend on the era considered, as the decay in 
the oscillations will depend on whether the tower is evolving during radiation or matter domination. 

Also, it must taken into account that the scalar fields that started oscillating during radiation domination 
will have a slight different behaviour during matter domination than the ones that start oscillating during 
matter domination. 

For radiation domination, pdm is simply given by: 

PDM = G. (^3/2-25^-372 _ ^-(3-25)j (jg) 

while for matter domination, two distinct cases must be taken into account depending on the value of a. 

If (7 < l/G, the evolution of pdm is described by: 

Pdm = J rn~'^^t~‘^dm = 

while for cr > l/tej it must be taken into account that some of the scalar fields in the CTSF have started 
oscillating during radiation domination, which means the integral to determine pdm must be split as: 


Pdm = G* X 


m ‘^dm = {m\ ^ 

+ 


( 12 ) 


5 













(13) 


Pdm = 


2C* cr3/2-26 

- 1 + 

3 — 46 \/mr 




(l- 26)(3- 46) 


l-2b,-2 _ _C^,-(3-2b) 
^ 1 - 26 


Pdm = 




3-46 


3/2-26 


/TOe 


6-2 - _^i-(3-26) 

1 - 26 / 1-26 


(14) 


with me = l/6e- 


In order to understand how well the previously obtained solutions describe the evolution of the CTSF, one 
can compare the results obtained with the integration of the analytic solutions presented in Q , considering 
the same set of scalar fields in both cases, which leads to figured 

p p 




Figure 2: In the left panel, the evolution of pT is displayed for 3 different sets of 6, the initial parameter, while the right panel 
shows the same quantities during matter domination. The values of b considered in the plots are the following, b = 0.2, b = 1.1 
and b = 1.45, with the decay of pT slowing as the value of b is raised. In both plots, the thick line denotes the approximation 
considered in this section, while the dashed lined refers to the analytic solution obtained in [^. 

If one aims to take this comparison further, it can be useful to compute how the equation of state, wt^ 
changes with the initial parameter b in both described cases, which leads to figure [3l This parameter is 
in fact a way to parameterise the growth or decay of pr, given that pr oc and therefore, by 

comparing the quantities for both cases, it is possible to understand how similarly px varies in both cases. 


w 


w 



Figure 3: In the left panel it is shown how the equation of state for the CTSF depends on the initial parameter 6, w{b), during 
radiation domination, while the right panel shows the same quantity during matter domination. In both plots, the thick line 
denotes the approximation considered in this section, while the dashed lined refers to the analytic solution obtained in [^. 
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Again, in both cases, only small differences are observed, meaning the simplified solutions are indeed a 
good match to the analytic solutions for the CTSF as well. 

It must be mentioned that the previous results for pT are only valid for tint = l/cr. 

If tini < I/cr, all the scalar fields in the CTSF will remain frozen until t = l/tr, and hence behaving like 
DE. This means pde = pdm = 0 for t < I/cr, considering the chosen set of scalar fields. 

On the other hand, if t^i > l/cr, all the scalar fields with masses a < m < rrii = IjUni will start 
oscillating at tini- This means there will be an extra amount of DM that must be taken into account, which 
will only affect the evolution of pdm- 

The new evolution for pdm can be obtained by replacing cr in the expressions for pdm, presented in 
eallll and 1141 with rrii and then adding the following term: 

-- ^ {£} ^ (£) <-) 
which corresponds to the energy of the extra amount of DM, taking into account it decays with a~^. 

Also, the results for w{h) presented in figure [3] will not depend on the set of scalar fields considered for 
the CTSF, given those are asymptotic results. 

However, extending the width of the distribution of the CTSF, in order to add an extra set of scalar 
fields with higher masses, can have a significant impact in the value of re(6) at a given time. 

In fact, if b is large, the only way the CTSF can account for DM as well as DE is to consider the previously 
mentioned extension, so that the DM sector will include more scalar fields. Given that extra set will have a 
significant contribution to pT, it is possible to obtain wih) 0 even for high values of 6, as mentioned in Q. 

If the CTSF is aimed to account for both DM and DE, the expression obtained for px in is actually 
better than the ones presented in this paper. 

As can be seen, the presented expressions for poM and pde both present singularities, depending on 
the value of b considered. The singularity presented in the expression for pdEi corresponding to 5 —>■ 3/2, 
is expected and was mentioned in the discussion on Q , but the singularities presented in the expression for 
Pdm, considering b 1/2 and b —>■ 3/4, are not physical and appear due to the approximation considered. 

However, the expressions determined in this paper are actually better when considering that the CTSF 
does not have to account for DM, given that it does not rely in the assumption that the DM sector of the 
CTSF dominates at tini- 

As a final note, despite TOs = 1/t in this paper rather than nis = H, as used in the cosmological 
results obtained previously in Q shall hold for this paper and the parameters shall be of the same order 
of magnitude with the present consideration for mg. Also, this new choice for TOs is solely due to it being 
better to simplify the calculations presented in this paper, but still a reasonable choice, as seen in the results 
presented in figure [1] 

3. Dark Matter - Dark Energy Interaction in a Continuous Tower of Scalar Fields 

As seen in the previous section, the CTSF can easily be split in two sectors, one that behaves like DE 
while the other behaves like DM. The mentioned splitting is possible due to the simplification considered, 
as it forces the scalar fields in the CTSF to be either in a DE or DM state, with the transition from a DE 
behaviour to a DM behaviour occurring instantaneously at t = 1/m. 

Since the scalar fields change from a DE state, when they are frozen, to a DM state, when they oscillate, 
the scalar fields are in fact transferred from the DE sector of the CTSF to the DM sector at t = 1/m. 
Therefore the transference of energy from one sector to the other can be seen as an exchange of scalar fields 
from the DE sector to the DM sector and hence the energy that is transferred at a given time will depend 
on the energy possessed by the scalar fields that change states. This is considerably different than the usual 
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interacting DM-DE models, in which the energy transferred at a given time depends in the energy of one or 
both of the whole fluids involved in the interaction. 

The energy transferred is given by the coupling Q, which will therefore depend in the energy of the 
scalar helds that start oscillating at a given moment. Said scalar fields will cover a small range in the mass 
spectrum, with m ~ TOs = 1/i, whose energy is given by eql?] 

Taking into account the normalisation due to the continuous distribution of scalar fields, Q can be written 
as: 


Q = (to)^ = ( 16 ) 

assuming that the CTSF is evolving during radiation or matter domination. 

Therefore, one finds Q oc during radiation domination, and Q oc during matter domi¬ 

nation, given that a (x as mentioned before. 

Switching back to time, the equations that take into account the interaction between the two components 
are given by: 


Pde = —Gift G 2b) ^ 

RDM + “iHpDM = , 

assuming that the scalar fields in the DE sector are frozen, which means wde 
that the energy flows from the DE sector into the DM sector. 


(17) 


— 1, and taking into account 


As can be seen clearly, both equations are not actually coupled, which means it is possible to obtain the 
evolution for rdm and pde independently. 

This is one of the biggest advantages of this kind of parameterisation, given that this is not the case 
with conventional models that consider Q to be parameterised by Q = SdeHpde + 5dmH poM, for which 
a more complex dynamical systems analysis must be considered in order to understand how pde and pdm 
evolve. 


Starting with the equation for DE, since it does not depend on the era considered, the following result 
is obtained: 


Pde 


f-(3-2b) 

3-2b 


(18) 


which fully agrees with eqE] obtained in the previous section. 

To obtain pDM^ the era in which the CTSF is evolving must be again taken into account, as the evolution 
of H will depend on the era considered. Therefore, for radiation domination rdm will evolve as: 


Pdm = Cit 


2C^, 

3-46 


while for matter domination the evolution will be given by: 


(19) 


Pdm = C2t 


j.-(3-2b) 

1 - 26 


( 20 ) 


The previous results obtained for rdm are in agreement with easfTT] and [HI obtained in the previous 
section, considering the simplified approximated evolution of the scalar fields, as long as: 


Cl 


2C* 3/2-2b 

3-46^ 

8 


( 21 ) 







and: 




C2 = 


3-Ab 


1-26 

3 / 2-26 


1-26 


m, 


1-26 


, a < l/te 
, a > l/te 


( 22 ) 


This means in fact the CTSF can be described as a model of interacting DM-DE during radiation and 
matter domination, as long as the adequate parameterisation for Q is considered. 


However, it must be pointed out that the effective value for wdm and wde can be quite far from the 
expected values (wdm = 0 and wde = —1), as can be seen in figure |H 


w 


w 



Figure 4: In the left panel, the equation of state for the two sectors of the CTSF, the DM and DE sectors, is plotted in 
function of the initial parameter 6, t6>(6), during radiation domination, while the right plot shows the same quantities during 
matter domination. In both plots, the thick line denotes the DE sector, while the dashed lined refers to the DM sector of the 
CTSF. 

Despite figure H] shows a wide range of possible values for wdm and wde, the scalar fields that compose 
each of the sectors will behave either as DM or DE at an individual level, which explains the name given to 
the different sectors. 

If one takes into account the equivalence to the CTSF shown and use the constraints on the initial 
conditions obtained in it is easy to understand which values are in fact reasonable for wdm and wde- 

As seen in for the CTSF to be in agreement with the observations, the value of h must be close to 1.5, 
which implies wde ~ —1, but also wdm ~ —1, if one considers a distribution for the CTSF with tini = Iju. 
Said assumption was used to determine the constants Ci and (72, presented in eqs. [H] and [H] respectively, 
but unfortunately does not allow for the whole system to account for the whole dark sector. 

If one wants for the whole system to account for the whole dark sector, an extra amount of DM is 
required, which will prevent the asymptotic behaviour presented in figure 0] to be reached until now and 
keep Wdm ~ 0. Said extra amount of DM can be obtained by raising the constants Ci and (72, which in 
fact corresponds to consider a larger width to the distribution in the CTSF. 

It must also be pointed out that, given Q does not depend on pdm, a change in Ci and C 2 will not affect 
the evolution of poM or pde, unlike what happens in other models of interacting DM-DE described in the 
first section. 

It is worth mentioning that this model does not allow for further constraints to be applied to the 
parameters on the CTSF model. 
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Unlike the models mentioned in the introduction, calculating the ratio R = pde/pdm does not make 
sense in this model, given that it breaks down when the DE sector of the CTSF starts dominate. This way, 
it is not possible to obtain constraints to the interaction, that would necessary lead to constraints on the 
parameters of the CTSF model. 

Therefore, it is also not possible to use this model to solve the ’’coincidence” problem due to the same 
reason, which was expected before hanih given that the asymptotic behaviour of the CTSF is equivalent 
to a cosmological constant, as found in [6|, which means an interacting model using the CTSF could never 
solve said problem. 

As a final note, it is relevant to mention that it would have not been possible to establish this equivalence 
using the expressions presented in Q for the different p’s, given that it was assumed that the whole energy 
of the CTSF is on the DM sector at tini . Despite this is a good approximation when the CTSF accounts for 
both DM and DE, this approximation prevents any energy from being transferred to the DM sector, which 
is crucial to establish this equivalence. 

4. Conclusions 

As aimed, this paper presents a new parameterisation for Q, the coupling between DM and DE, physically 
supported by an equivalence between the CTSF, first introduced in and a system of interacting DM-DE. 

This analogy is only possible if the scalar fields belonging to the CTSF are assumed to change from a 
DE behaviour to a DM behaviour instantaneously. As seen before, this is only reasonable for radiation and 
matter domination, but not when the DE sector of the CTSF is the dominant component of the universe. 
Despite this drawback, the results obtained for the CTSF considering this approximation for the scalar fields 
pretty much match the results obtained using the analytic solutions for the scalar fields presented in 

By splitting the CTSF in a DM and a DE section and treat it as a DM-DE interacting system, it is 
possible to obtain the same results as before, as long as it is considered that the DM fluid has w = 0, while 
w = — 1 for the DE fluid, and also that Q corresponds to the energy of the scalar fields that start oscillating 
at a given time and hence change from a DE state into a DM state. 

This way it is clear the CTSF is basically equivalent to a system of interacting DM-DE during radiation 
and matter domination, which introduces a new perspective to how the CTSF actually works. 

Furthermore, this equivalence allows the introduction of a different parameterisation for Q that does not 
depend on any properties of the DM or DE fluids, but rather on the energy of the individual scalar fields 
that constitute the DM and DE fluids. The energy is therefore quantised in the scalar fields that constitute 
the CTSF, giving a whole new perspective on how the interactions between DM and DE can work. 

Despite this paper focus only on the CTSF model and how it can be seen as an interacting DM-DE 
model, it is possible to generalise the calculations done to fit other models involving multiple scalar fields. 

Once again, it is important that the transitions a scalar field undergoes from one state to another occur 
quickly, so that it is possible to approximate said transition to be instantaneous if the interacting fluid 
approach is to be applied to other systems with multiple scalar fields. However, it is not necessary that the 
scalar field can only be in a DM or DE states for this formalism to be applied, as long as it is taken into 
account the actual value of w for each of the state the scalar field can be at and that the expression for Q 
is modified so that it corresponds to the actual energy a scalar field carries when switching states. Given 
that Q will not depend on the total energy of DM or DE, this also means it is much easier to understand 
the evolution of poM and poE without resorting to a dynamical systems analysis. 

It is also relevant to point out that the mass in the multiple scalar helds system does not have to be 
continuous. However, a discrete distribution over the mass yields discrete transferences in the energy from 
one component to the other. This would be relevant for the model presented in Q, for instance. 

Last but not least, it is also relevant to mention that the approximations introduced for the evolution of 
can be quite useful, given that it allows for a simple and accurate expression to be obtained for pT, as 
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seen in section 2. Despite the expressions obtained end up not being as simple as the ones obtained in 
these expressions are valid when the DM sector of the CTSF considered is not clearly dominant a.t t = to, 
which is not the case with the expressions presented in Q, given those rely on pdm ~ Pt at the initial time. 
Therefore, any future project featuring the CTSF can benefit from this simplification as it will speed up 
significantly any computational calculations that are required for radiation and matter domination, given 
that it is possible to use a simple analytic expression rather than a considering the numerical integration of 
analytic solutions for (j)m- 
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